function foldPlot(dataFileName, dataFilePath, title)
if (nargin < 1)
    disp(['need at least 1 arg']);
end
if ( nargin < 2)
    dataFilePath = './';
end
if (nargin < 3)
    title = 'Title';
end

%clc;
disp(['']);
disp(['====================================================']);
disp(['= file name = ',dataFileName]);
disp(['= file path = ',dataFilePath]);

dataFile = strcat(dataFilePath,dataFileName);
disp(['= datafile = ',dataFile]);
% import data file
foldData = importfile(dataFile);

nbin    = foldData.data(:,1);
phase   = foldData.data(:,2);
mean    = foldData.data(:,3);
stdev   = foldData.data(:,4);
nPoints = foldData.data(:,5);


% fit the data
[coeff, fitStats] = cosFit(phase, mean, stdev, nPoints, dataFileName, title);
%disp(coeff);
%disp(fitStats);
%save coeff;
%save fitStats;

% sse = sum of squares due to error
% rsquare = corefficient of determination
% dfe     = defrees of freedom
% adjrsquare = degree of freedom adjusted for coefficient of determination
% rmse = root mean squared error (standard error)

